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The one-dimensional isotropic quantum Heisenberg spin systems with random couphngs and 
random spin sizes are investigated using a real-space renormalization group scheme. It is demon- 
\ strated that these systems belong to a universality class of disordered spin systems, characterized 

by weakly coupled large effective spins. In this large-spin phase the uniform magnetic susceptibility 
diverges as with a non-universal Curie constant at low temperatures T, while the specific heat 
O . vanishes as T''| InT] for T ^ 0. For broad range of initial distributions of couplings and spin sizes 

■ the distribution functions approach a single fixed-point form, where 5 « 0.44. For some singular 

initial distributions, however, fixed-point distributions have non-universal values of 5, suggesting 
, that there is a line of fixed points. 

(N 
o 

I. INTRODUCTION 

. Over many decades one-dimensional (ID) quantum spin systems ('quantum spin chains') have attracted a lot of 
c/3 ' interest and led to the deveJopmcnt of many theoretical methods which are now commonly used for the study of other 
^ highly correlated systems.El Despite the apparent simplicity of quantum spin chains, they show a wealth of physical 
properties which give a key to om: understanding of various phenomena, e.g., quantum phase transitions, topological 
^ order, and fractional statistics. Elij Since the discovery of various quasi-lD materials, the study of ID spin systems, 
which is mainly based on the Heisenberg naodel and its variatiwis, is also of experimental relevance. Examples of such 
materials include so-called NINO, NENP,Qtl and SrsCuPtOgfl In particular, the latter system belongs to a class of 
Q i compounds which is compositionally very flexible and has been under intense experimental investigation over the last 
^ ' few years. This type of quasi-lD system was first reported in Sr4Pt06 by Randall and Katz,Ll and it is now possible 
to produce compounds of the form SrsMNOe in various combinations with M = Cu, Mg, Zn, Yb, Na, Ca, Co and N 
^ = Pt, Ir, Rh, Bi. 

^ , Disorder effects play a particularly important rolR-|in ID quantum spin systems, as even small deviations from 
\^ • the regular system often destabilize the pure phasesJj Real experimental systems naturally contain impurities and 
' other types of disorder. Therefore it is very important to understand the influence of disorder on the properties of 
. such systems in order to interpret experimental results. To our knowledge, the first ID spin Sjystem recognized for 
its disorder belongs to the class of charge-transfer salts TCNQ (tetracyanoquinodimethanide) .□ These systems have 
been successfully described by a ID spin-i Heisenberg model with random strength of antiferrornagnetic exchange 
couplings between the spins. A more recent example of disordered spin chains is Sr3CuPti_2,Ir2;06.li3 While the pure 
compounds SrsCuPtOg {x = 0) and SrsCuIrOe {x ~ 1) are antiferromagnetic (AF) and ferromagnetic (FM) spin 
systems, respectively, the alloy SrsCuPti-ajIrajOg contains both AF and FM couplings. The fraction of FM bonds is 
^ , simply related to x, the concentration of Ir ions. In a previous work, we modeled SrsCuPti-jjIrjjOg with a nearest- 
I ' neighbor spin-i Heisenberg chain, where the exchange coupling between neighboring spins is + J or — J with probability 
' O ' p and 1 — p respectively. The methods we used (high-temperature expansion and transfer matrix approximation) give 
^ reliable results down to J/5, where T is the temperature and ks is Boltzmann's constant. In this regime 

the numerically calculated magnetic susceptibility x{T) and the specific heat C(T) interpolate smoothly between 
the corresponding nuantities of purely FM {p = 1) and AF {p = 0) chains, giving good qualitative agreement with 
experimental data.lliJIlj At temperatures below ~ J/ks the effects of disorder become significant. We demonstrated 
that in this temperature regime spins correlate within AF and FM segments of the chain separately. The emerging new 
5J] ' degrees of freedom which dominate the thermodynamics are (large) effective spins each corresponding to a correlated 
[ segment. The size of these spins and their residual interaction are set by the local disorder and hence is random. From 
exact diagonalization of finite segments we concluded that the low-temperature physics of the random spin system is 
described by the effective Hamiltonian, 

'H — ''^JiSi ■ Si-f-i , (1) 

i 

where both the couplings Ji, which may have either sign, and the spin sizes Si are random. In particular, we emphasize 
that the resulting distribution of Jj in Eq. (0) is broad and dense, in contrast to the discrete distribution of the initial 
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model. In this paper we take Eq. (|^) as starting point. We discuss the low-temperature propierties of this model, and 
the various fixed points encountered for different initial distributions of couplings and spins.tJ 

Before going into details we briefly summarize our results and the method we use, which is a generalization of the 
RSRG scheme introduced by Ma, Dasgupta and Hu (MDH) in 1979 to study the ID spin-i random antiferromagnet 
(RAF). til The RAF, where all couplings ^e antiferromagnetic but vary in magnitude, has also Jaeen investigated 
using Kadanoff block spin RG techniques,E3 and more recently by the density matrix RG methQd.EJ The method of 
Dasgupta et al. has proven to be the most successful one, and was recently extended by Fisher Jl3 who solved the RG 
equations exactly. 

In the MDH RSRG scheme, a decimation of degrees of freedom occurs through the successive formation of spin 
singlets from the most strongly coupled spin pairs. This scheme conserves the form of the Hamiltonian in the original 
model, but changes the distribution of couplings, which gradually approaches a flxed-point form. The model in 
Eq. (|l|) contains arbitrary spin sizes and couplings with random sign, so that, in general, two correlated spins do 
not combine into a singlet. Rather they form an effective spin with renormalized couplings to its neighbors. Here 
we introduce a modified RSRG scheme which takes this into account. Like the MDH RSRG scheme it conserves the 
form of the Hamiltonian, Eq. (|l|), but changes the distributions of couplings and gaps. The RG fiow generated can 
therefore be thought of as a flow in the space of distributions of couplings and gaps.EJ We demonstrate that for a wide 
range of initial distributions of couplings and spins, the RG flow of the distribution functions eventually approaches 
a single universal fixed point. This fixed point represents the following physical properties. Both entropy and specific 
heat vanish as InTj, where the power a ~ 0.22 is rather small. The exponent a appears also in the non-linear 
magnetization where M{H) cx for sufficiently large fields H. For very singular initial distribution of the 

couplings we find that a takes a non-universal value, suggesting the presence of a fixed line. We also find a surprising 
fact that, for both the universal and non-universal fixed points, the susceptibility follows the Curie behavior down to 
zero temperature with a non-universal Curie constant. We finally shoffi-Jtiaw the two previously known random phases, 
the random singlet phase (RSP) and the random dimer solid (RDS)E3'c2l, both of which correspond to Eq. (|l|) with 
all Si = 1/2 and all Ji > 0, are unstable against the admixture of an arbitrarily small concentration of FM couplings 
and/or larger spins. 



Our paper has the following structure. We start with a brief review of t he MDH RSRG scheme in Sec. [I A, before 



we generalize it to chains with both AF and FM couplings in Sec. |IB|. In Sec. |I^ we analyze the distributions 
of spins and couplings, and their scaling forms close to a fixed point. Wc perform the RG scheme numerically 
by simulating random spin chains with various initial distributions of couplings and spins. The numerical results 



shown in Sec. IV confirm the scaling forms conjectured in Sec. [II, but also reveal that random spin chains with 
very singular initial distributions of gaps fiow to non-universal fixed points. In Sec. ^ we derive the scaling forms of 
thermodynamic quantities, and in Sec. [v| we comment on the approximations involved in the RG transformation. 



Finally we summarize our results in Sec. VII, and compare the large-spin phase (LSP) to other disordered phases. 



We also discuss the stability of the various phases and in particular the RG flow between the LSP, RSP, and RDS. 

II. THE RENORMALIZATION-GROUP SCHEME 

To study the low-temperature pmperties of systems described by Eq. (P, we generalize a RSRG method introduced 
by Ma, Dasgupta and Hu (MDH).Efl We start with a brief review of the MDH scheme for the RAF with Si = 1/2 and 
random Ji > 0. 

A. The MDH RG for antiferromagnetic spin-i chains 

Consider an antiferromagnetic nearest-neighbor Heisenberg spin-^ chain in which the largest coupling is Jq and the 
remaining couphngs Ji are distributed according to P(Jo; Ji)- We focus on the link with the largest couphng, Ji — Jq, 
and the terms in the Hamiltonian (Q) that involve the spins and S^+i, (see Fig. ^), 

n' ^H'o + n'l (2) 

where 

Hq = Jo^i ■ Si+i , 

Ti'i = Ji-iSi_i • Si + Ji+iSi+i • Si+2 . 
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If the distribution P(Jo; J) is broad, Ji±i are typically much smaller than Jq and we can treat Ti'j as a perturbation 
to H'q. In the ground state of H'q the spins Si and Si+i form a singlet, and the energy gap to the excited states is 
Jq. This ground state is four- fold degenerate since the unperturbed Hamiltonian H'q does not involve the directions 
of Si-i and Si+2- H'j lifts the degeneracy and splits the unperturbed ground state into a singlet and a triplet, and 
the low energy spectrum of the four-spin Hamiltonian (|^) is described by an effective Hamiltonian 

Ti^^ = JSj_i • • (3) 

The effective coupling J is determined from the energy splitting of the unperturbed ground state, and to second order 
in Ji±i/jQ the result is J/Jq — '''Z^ji'*'^ . Physically the weak interaction between Si_i and 8^+2 is mediated by 

exciting virtual triplet states in the interjacent spin pair. 

In the Hamiltonian (]l|) we replace the terms in H' with the effective interaction H"^ in Eq. (||) to get an effective 
Hamiltonian for the low-energy degrees of freedom of the spin chain. Repeating this procedure and successively 
replacing the strongest remaining coupling in the chain preserves the form of the Hamiltonian but changes the 
distribution of couplings and, in particular, lowers Jq, the largest remaining coupling in the chain. If P(Jg, J) 
is the distribution of couplings at a point when the lacgest remaining coupling is Jq, then the removal of bonds 
Ji G [Jo — dJo, Jo] generates a flow equation for P( Jo, J)lij 

'^'^^"''^^ = -P(Jo; Jo) / dJiJJ2 PlJo; Ji)P(Jo; J2)5(J - Ji J2/2 Jo) . (4) 
"-'0 Jo 

The flow equation is derived under the assumption that there are no spatial correlations among the bond strengths. 
This is indeed the case if there are no correlations in the distribution of couplings in the initial chain. It has been 
shown that if the initial distribution of bonds is normalizable, Eq. (jj) has,|a,unique fixed-point solution that governs 
the low-energy physics of random bond antiferromagnetic spin-i chains .EZrEH 



B. Generalization of the MDH RG 



We apply the same strategy to the random spin chain with couplings of either sign and random spin sizes. In 
contrast to the previous case, a link is determined not only by the coupling strength but also by its left and right 
spin, {Ai, S'i, S'i_|_i} (see Fig. ^3). We define as the energy gap between the ground state multiplet and the first 
excited multiplet in the corresponding two-spin Hamiltonian Ti, — JiSi ■ S^+iH 



I |J.|(S.-f S.+i) 



M\s.. 



1) 



Ji < (ferromagnetic link) , 
Ji > (antiferromagnetic link) 



(5) 



We assume a broad distribution of interaction energies and focus on the link in the chain which, if completely isolated 
from the rest of the chain, requires the largest energy A = Ao to excite the ground state multiplet. We consider 
the situation illustrated in Fig. |l]c, where {Aq, Sl, Sb} is the strongest link in the chain, and {Ai, Si, Sl} and 
{A2, Sr, S2} are its adjacent Hnks. In the spirit of the MDH RSRG scheme we replace the strongest link {Ao, Sl, Sr} 
with an effective spin of size S = {Sl^SrI representing the lowest-energy multiplet of maximum (minimum) spin for a 
ferromagnetic (antiferromagnetic) link. The residual effective interaction for Si, S, and S2 is calculated perturbatively 
in £1,2 = ^1,2/^0- The effective interaction is isotropic, and to first order in £12 given by 
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where 



fla a a \ — Sl{Si+Sr — Sl) 



{Sl — Sr 


+ l){Si + 
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U{Si,Sl, Sr) - — — — — -r^ (8d 

[dr — 5l + + ol) 

„ p c \ Sl{\Si - Sl - Sr\ + 1) . 
MS,, Sl, Sr) - + + (8e. 

J- /r. c C \ Sl{Si+Sl+Sr) 

{SL + SR)iS, + SL) ■ 



A derivation of these equations is shown in Appendix Al. From the knowledge of the gap Ai and the sign of Ji, 
Ji is readily calculated via Eq. (||). Similarly A2 is obtained by replacing 5*1 by S2 and Sl by Sr in Eqs. @) and 
(^. These equations do not require the spins to be multiples of 1/2, and from now on we regard spins as continuous 
variables. The case where the strongest link is antiferromagnetic with 5*^ = Sr is not accounted for in Eq. (^. In 
this case the two spins Sl a nd 5 *^ form a singlet, and the leading order contribution to the effective coupling between 
Sl and 5*2 is (c.f. Appendix |A 2| ) 

J-^SLiSL + l) , (9) 

which is easily translated into A via Eq. (||). For Sl = Sr = 1/2 this is the original MDH RG transformation. 

As in the original RSRG scheme, the effect of successively forming effective spins is to change the distributions of 
gaps and spins (links) without changing the form of the Hamiltonian. In analogy to the probability distribution for 
the couplings in the RAF, we define probability distributions of ferromagnetic (P^) and antiferromagnetic (P^) links 
where the largest remaining gap in the chain is Aq, 

P^{Ao;A,Sl,Sr) , (10a) 
P^iAo;A,SL,SR) . (10b) 



The probability distributions P^^'^ are symmetric in Sl and Sr, and obey the normalization condition 

pAa poo 

/ dA / dSLdSR [P^(Ao; A, Sl,Sr) + P^iAo; A, Sl, Sr)] = 1 , 
Jo Jo 



(11) 



for any value of Aq. From P^ and P^ we can calculate the distributions of spins, gaps and coupling constants. As a 
special case, the original spin-i antiferromagnetic chain studied in Ref. |l^ corresponds to 

P^{Ao;A,Sl,Sr)^S{Sl~^)S{Sr~^)P{Ao;A) , (12a) 
P^ {Ao; A, Sl,Sr)=0 . (12b) 

If there are no correlations between neighboring links (except for the obvious correlation that they share one spin), 
the flow equations for P^^^ are 

^^Fi[P^,P^] , (13a) 
dAo 

'^-F,iP^P-i . (m) 

which generalize the MDH RSRG flow equation (^. In Eqs. (^3|) Fi and F2 are two (non- linear) functionals of P"^ 
and P^ , whose explicit forms depend on the functions /„ in Eq. (ph, c.f. Appendix H. 
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III. SCALING FORMS OF THE PROBABILITY DISTRIBUTIONS 



As links are replaced by effective spins, the effective couplings and, hence, the gaps of the links decrease. At the 
same time the average distance nao between neighboring effective spins as well as the magnitude of the effective spins 
increase. (Here oq is the original lattice constant and n is the ratio of the number of original spins to the number of 
effective spins.) We expect that the link distributions eventually approach a fixed point where P^, and n exhibit 
scaling behavior 

P^(Ao;A,^.,^«) = Ao--g-^(^,^,^) , (14a) 
P^(Ao;A,5i,5«) = A-^-Q^f-^,-A^,-^^ , (14b) 



and 



n^A-^ . (15) 



The exponent 5 is related to the dynamical exponent z {z — 1/5). The seven exponents in Eqs. ( p^ ) and (|T^) are not 
all independent. Indeed, we argue below that 

ap = OLA = a (16a) 
/3f = /3a = 1 (16b) 
7F = 7A = 1 — 2q! (16c) 
5^ 2a (16d) 



so that in the scaling regime 

P^^ 

A, 



- Ti^t?^(^/^o, ^lAo", ^flA^) (17a) 







= Ti^Q''(A/Ao, 5lA;?, SrIS.'^) (17b) 

^0 



with length scaling as 



Ao'" . (18) 



The relations (Eq) have been confirmed in numerical simulations (c./. Sec. IV) and can be understood as follows. 

Let X = Na/JNa + Np) be the fraction of AF links in the chain. Both FM (x = 0) and AF {x = 1) chains are 
unstable towards a small concentration of couplings of the opposite sign. To see that a; = 1 is unstable, we note 
that unless the effective spin formed is a singlet, the removal of a link in an AF chain converts one neighboring link 
into a FM link. Similarly a; = is unstable because an isolated AF link in a FM environment always survives (the 
removal of a FM link does not change the signs of its neighboring links, and if the AF link itself is removed, one 
of its FM neighbors is converted into an AF link). This implies that for small enough x the absolute number of 
AF links is constant, so that the fraction x of AF links increases as links are removed. Thus, unless we start with 
a completely FM random spin chain or a purely AF random spin chain with uniform magnitude of the spins, the 
fixed-point distribution contains both FM and AF links. Having established that both xq and 1 — xq, t he fraction 



of AF and FM links at the fixed point respectively, are non-zero, we can easily derive the relations (16a-c). Since a 
finite fraction of the spins belong to both a FM and an AF link, there cannot be a separation in scales between spin 
sizes in AF and FM links, i.e., aA — ctp — a- (3f — ^ {Pa — 1) follows from the fact that the average gap in the FM 
(AF) distribution, when measured in units of Aq is finite and independent of Aq. 7a = 1 — 2a follows trivially from 
the finiteness of xq 



Xo 



i-Ao poo 

/ dA / dSLdSRA-^^Q^{A/Ao,A^SL,A^SH) 
Jo Jo 

^l-2a-7A / / ds'ds"Q'^{A',s',s") , (19) 

Jo Jo 
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which must be independent of Aq. An analogous argument for 1 — xq impHes 7^? = 1 — 2a. 

The last relation (16d) requires a more detailed analysis of the way effective spins are formed by correlating the 
original spins in clusters. The size of the effective spin corresponds to the spin quantum number of the ground state 
of the cluster. Since the spin system is not frustrated the spin quantum number is determined from the classical 
correlation of the spins (parallel and antiparallel). The total spin of a cluster of n spins is then given by 



S 



(20) 



where two neighboring spins enter the sum with the same (opposite) sign if their mutual coupling is ferromagnetic 
(antifcrromagnetic). This leads to a typical random walk problem which results in the scaling 



,1/2 



(21) 



From this we conclude 5 = 2a. 

The remaining independent exponent, which we take to be a, is related to the average renormalization of the gaps. 
Consider two neighboring effective spins with a coupling corresponding to a gap A. Suppose also that each of the 
two effective spins is made up of 2*^ spins at some larger energy scale A'. As the energy scale is lowered from A' to 
A and the 2*^ spins form one large effective spin, the gap A is typically renormalized 2k times (fc times in the course 
of formation of the left effective spin and similarly k times from the right). If the magnitude of a gap is reduced on 
average by a factor a each time a neighboring link is replaced with an effective spin, then 



n'2'' =n' (a2'=)~ 



from which we read off the relation 



In 2 

4 In a 



(22) 



(23) 



Using the scaling form, we can get some information on long-distance behavior of spin-spin correlation functions. 
Let us introduce two kinds of spin-spin correlation functions which characterizes the correlation between spins at low 
temperatures. The first one is the usual spin-spin correlation function, 



Ci(*-j) = (S.-S,), 



(24) 



where ( ) represents both thermal average and average over random configurations. Since the number of AF bonds 
between Si and Sj is random, the two spins may either be in parallel or antiparallel. Thus, after taking the random 
average the correlation function decays exponentially for large |« — j| even at zero temperature: 



Ci{i-j) (xexp< 



In 



1 



|2p-l| 



+ iTre{l - 2p) 



(25) 



where p is the density of FM bonds. Therefore this correlation function does not reflect the correlations leading to 
the formation of effective large spins. An appropriate correlation function is 



(26) 



where r/ij — Yit—i ^S^^{~Jk) for j > i. At finite temperature this correlation hmction should also decay exponentially 
for large |j — J| with the correlation length n: 



C2{i-j) cx exp(-|i - j|/n) 



for T > 0. 



(27) 



From Eq. (|l8|) we find that the correlation length grows with decreasing temperature as n cx T It is likely that 
at zero temperature the correlation function decays algebraically as 



C2{i- j) oc 



1 



(28) 



Unfortunately we cannot determine the value of this new exponent v from our numerical RG scheme. 
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IV. NUMERICAL RESULTS 



We perform our RG scheme by numerical simulations. We start each simulation by generating a chain according 
to independent probability distributions for gaps and spins. In each decimation step we pick up the strongest link in 
the chain, replace it with an appropriate effective spin, and renormalize the neighboring bonds. To keep the number 
of links fixed, one site is finally added in one end of the chain. This procedure is then iterated until the shape of the 
distribution function of links no longer changes. In this way we have iterated sixteen chains with both non-singular 
and singular initial distribution, see Tab. I 

In all our simulations the distributions of links eventually converged to some fixed-point distributions. The distribu- 
tions rather quickly take the rough forms of the fixed-point distributions, while the final approach and, in particular, 
the convergence of the exponents in Eqs. (^6|) to their final values are very slow and take place over up to five orders 
of magnitudes in length (ten orders of magnitude in energy). Our numerical simulations demonstrate that unless the 
initial distribution has a high degree of singularity for small gaps, the distribution of links in the chain eventually 
flows to a universal fixed-point distribution of AF and FM links. If the initial distribution of gaps is more singular 
than P{A) ~ A"^'', yc ~ 0.7, our numerical simulations suggest that the corresponding fixed-point distribution is 
non-universal. This is analogous to the RAF where it has been shown that extremely singular components in the gap 
distribution are conserved in the RG flow.tO In the case of the RAF yc — 1, so that the condition for a chain to flow to 
the universal flxed-point distribution coincides with normalizability. In contrast to the RAF, |tiiere may be p hysical 
situations where a random AF/FM spin chain flows to a non-universal flxed-point distribution,E3 c.f. Sec. |Vl4 Below 
we summarize the numerical results in the case of regular and singular gap-distributions. 



A. Regular and weakly singular distributions 

If the initial distribution of gaps is regular or at least less singular than ^'(A) ^ A^'^'^, we flnd that the link- 
distributions in all chains we have studied (chains A, B, C, D, F, H, and I in Tab. |) eventually converge to the same 
universal distribution with the characteristic scaling form ( p^ ) proposed in section II] . The flxed-point distribution 
functions arc illustrated by various cross-sections in Figs. || and |[ We flnd that the ratio of AF links stabilizes around 
X = 0.63, thus confirming the conjecture in section III that both xq and I — xq are non-zero. The exponents aA,F and 
/3a, F in Eqs. ( p^ are deduced from the scaling of the averages (S) and (A) with Aq in the scaling regime. This is 
illustrated in Fig. ^ for chain C in Tab. |, where the average gap and spin size are plotted versus Aq in a log-log plot. 
Similarly the ratio n of the number of original spins to the number of effective spins is plotted versus the maximum 
gap Aq in Fig. and the evolution of the exponents as functions of Aq are plotted in Figs. ^ and ^b. We find the 
fixed-point values of the exponents to be 

aA,F = 0.22 ±0.01 , 
(3a,f = 1.00 ±0.005 , 
(5 = 0.44 ±0.02 . 

The exponents in the FM distribution agree with the ones in the AF distribution within numerical accuracy. Thus 
our numerical findings confirm the scaling forms ( p^ ) and ( [l^ ) for the probability distributions of the links as well as 
the relations between the exponents, Eqs. (p^. Identifying a with either ua, cif or gives consistently 

a = 0.22 ±0.01 . (29) 

An interesting observation we have made is that the ratio 2a/ 5 generally stabilizes to its fixed-point value of 1 before 
the two exponents S and a separately converge to their corresponding fixe d-po int values. This confirms the robustness 



of the 'random walk' argument in Sec. [II leading to the relation in Eq. ( 16d ). 

The typical expansion parameters in the perturbative calculation of the renormalized gaps are the median ratios 
between A"^'^ and Aq. At the fixed point these ratios are 

A^/Ao^O.2 , 
Af/Ao^O.3 , 

where we denote the median gap in the distributions of antiferromagnetic and ferromagnetic links by A^ and Af 
respectively. As expected from the increase in effective spin size, the formation of a singlet on a link {Sr = Sl) 
becomes increasingly rare as the fixed point is approached. 
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B. Strongly singular distributions 



For chains with a very singular initial distribution of gaps the convergence to some fixed-point distribution is 
generally even slower than for regular (or modestly singular) chains. Furthermore, in these cases the fixed-point 
distribution that is eventually approached as well as the scaling exponent appears to be non-universal. This is 
illustrated in Fig. ^ where the fixed-point distributions for four chains are plotted. The dotted distributions correspond 
to chains where the initial gap-distributions are P(A) c>c with x = |3/4 , 4/5, 7/8} (Tab. |). The solid curve is 
the corresponding fixed-point distribution for regular chains discussed in Sec. |IV A . For reasons discussed below, the 



actual form of the distributions corresponding to singular initial distributions may be quantitatively incorrect, but 
they deviate clearly from the fixed-point distribution of regular chains. Numerically we find that chains with initial 
gap-distribution more singular than 

P(A) - A-y- with 0.65 < 2/c < 0.75 (30) 

flow to non-universal fixed-point distributions. From log-log plots of the universal fixed point-distribution of gaps. 
Fig. H, we find that the distribution of FM gaps diverges as P^(A) ~ /\-o.44 ^^yq distribution of AF gaps as 
P^(A) ~ /\~o.7o gjj^g^u gaps. Thus, even from a regular link-distribution the RG transformation itself produces 
a singular fixed-point distribution of gaps, where the degree of the singularity, P(A) ~ A^^'=, is set by details 
in the RG transformation rather than by the initial conditions. If a distribution is more singular than A~^% the 
singular component is conserved in the RG flow. This is supported by our numerical results that yc in Eq. (^0|) is 
somewhere between 0.65 and 0.75, which agrees with the singularity in the universal flxed-point distribution of gaps 
P(A) ~ /\-o.70_ -yyg conjecture that in chains with gap-distributions more singular than A~*''=, the low-energy fixed 
point is not determined by the RG transformation alone, but also by the singular distribution of extremely weak links. 
In these chains we expect the fixed-point distribution as well as the value of the scaling exponent to be non-universal. 

The picture we present is in close anabgy to the RAF where Fisher has shown that the flow equation (||) conserves 
very singular components of P(Ao, A).EZl In the case of the RAF y^ = \, implying that any normalizable (and hence 
physical) initial distribution eventually flows to a universal fixed-point distribution. 

The fact that a very singular fixed-point distribution is dominated by the weakest links in the initial chain casts 
some doubt on the numerical results for such chains at low energies. Indeed, since all our chains have only a finite 
number of links, the singularity cannot be resolved perfectly. Hence the number of initially extremely weak links 
which are important for the low-energy behavior of the chain are relatively few and we do not expect the numerical 
results for the extremely singular chains to be quantitatively correct. This explains why the singularity seems to 
soften at very low energies while we argue that it should remain constant. Still, the qualitative result that regular (or 
slightly singular) chains flow to a universal flxed-point distribution while more singular distributions do not, should 
be correct. 



V. THERMODYNAMICS 

A. Entropy and specific heat 

The scaling forms in Eqs. ( p^ ) and ( p^ ) allow us to determine the universal temperature dependence of various 
thermodynamic quantities which may be measured in experiments. Let us start with the entropy and the specific 
heat. At finite temperature T the renormalization group fiow stops at Aq ^ ksT due to thermal fluctuations which 
prevent the formation of even larger effective spins. At this point, all pairs of spins in links with gaps larger than 
Aq ~ fcsT form large effective spins. Since the distribution of gaps is broad, the interaction energies between the 
effective spins are typically much smaller than Aq kgT, and each large spin moves essentially independently. The 
entropy per unit length is hence 

a(T,^ = 0)^fc.ln(2(5e.) + l)^^,.,^^^, (31) 
L n 

in the scaling regime where {Scs) ^ 1- Note that the assumption of independent effective spins leads to an overestimate 

da 
dT 



of the entropy. From the relation C{T) = follows the specific heat per unit length 



C{T,H = 0) 



cxr^"|lnr| . (32) 
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This is qualitatively different than in the random spin-i antiferromagnet where a{T) cx |lnT| ^. This is also in 

contrast to the uniform ID antiferromagnet where aApiT) (x T and the uniform ferromagnet, crpM (T) cx VT. The 
fact that both the entropy and the specific heat of the random-exchange spin chains go to zero with a rather small 
power reflects the presence of large number of uncorrelated spin degrees of freedom at low temperature. 



B. Static magnetic susceptibility 

By analogous arguments, the essentially uncorrelated large effective spins give a Curie-like contribution to the 
magnetic susceptibility per unit length 



L SknT n 



A, 



-2a 



T A"2a 



c 
T 



(33) 



The T^^ Curie behavior is usually a signature of uncorrelated spins. We emphasize that this is not the case for 
the random spin chain. Rather, most of the original spins are strongly correlated, and the Curie-like temperature 
dependence follows from the scaling relation n ^ (S^ff), Eq. (]2l|) . The Curie constant c can be calculated in terms of 
the original spin distributions as follows. As discussed in Sec. \ni[ the magnitude of the effective spin representing a 
segment of n frozen spins is given by the sum 



q2 



(34) 



where the staggering factor Si is defined by 

Si+i = -6tSgn{ Ji) ; Si ^ I , 
and S'i > is the spin size of the elementary spin at site i. Averaging over the (initial) disorder, we obtain 



(35) 



\i=l / i^j 

Defining p as the probability for a bond to be ferromagnetic and using (SiSj) — {2p — 1)1*^^1 we get 

" 2p-l 



{St: 



1-p 



Oil 



and in the limit of small T (large n) we find the Curie constant to be 



3fcB 



{St) + 



2p-l 
1-P 



{Si^ 



(36) 



(37) 



(38) 



This result coincides with the low-temperature susceptibility we obtain for the analogous classical spin chain.0 Note 
that this low-temperature Curie constant is in general different from the high-temperature value 



It follows that c = c only if 



1-P 
2p- 1 



(39) 



(40) 



Thus, in the random spin chain we expect the magnetic susceptibility to cross over from one Curie-like behavior at 
high temperature (T > J/ks) to a different Curie-like regime at low-temperature (T < J/ks)- 
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C. Magnetization at finite H 



In a finite magnetic field H and at finite temperature T the RG flow is interrupted either by the thermal energy kgT 
or by the magnetic Zeeman energy Ezm = fJ'iSes)!!. If ksT > Ezm the chain is dominated by thermal fluctuations 
and the magnetization is given by xH ■ If ksT < Ezm the magnetic field drives the system away from the fixed 
point of zero magnetic field into a state of aligned effective spins where the magnetization eventually saturates. In 
this regime a non-zero magnetic field starts to align the effective spins at an energy scale Aq ~ /Lt(S'eff)ff. With 
above scaling properties this means Aq ~ H^^^^", so that the saturated magnetization per unit length becomes 
M/L ~ /i(S'cff)/7i ~ jja/i+a^ rpj^j, condition that the chain is not yet dominated by thermal fluctuations is ksT < 
Aq ~ Summarizing these arguments, we get 

M{T,H) (htt^ ■.T'+'-^bH 

L \ H/T : ri+" » bH, ^ ' 

where & is a dimensionful non-universal constant. Similarly the entropy goes rapidly to zero at T^+°' « bH when the 
magnetic fleld starts to align the spins. 



VI. COMMENTS ON THE RSRG SCHEME 



In this section we discuss the validity of our RSRG treatment and various approximations we used. As we have 
seen in the previous sections, the formation of an effective spin yields new interactions among the remaining spins. 
These interactions were calculated perturbatively, where the (average) perturbation parameter is e = At/Ao. To first 
order in e, only nearest neighbor Heisenberg terms are induced, and the functional form of the Hamiltonian is 
preserved in the RSRG transformation. However, the terms in (|l]) are not the only ones allowed by the symmetry, and 
in general we expect more complicated isotropic interactions to appear ifjiigher order corrections in e are included. 
In the original MDH RSRG, Aj/Aq ^ as the fixed point is approached, tD and the perturbative treatment becomes 
exact. In our case e stabilizes at a finite value around 0.2 at the fixed point. Thus we have to analyze here to what 
extent higher order terms can change our results. 

The basic assumptions of the RSRG scheme are that the two spins which are most strongly coupled to each other 
form one effective spin, and that any breaking up of this spin pair involves such a large energy that we can regard 
the effective spin as a rigid object. There are two criteria for these assumptions to be valid. First, the energy cost 
for breaking up the strongest spin pair, i.e., the energy gap A, must be much larger than the energy available in 
neighboring spin pairs. Second, non-nearest neighbor couplings have to fall off sufficiently rapidly with distance, so 
that many weak couplings cannot accumulate sufficient strength to break up the spin pair. The second criterion is 
essentially equivalent to the absence of strong frustration in the system. Below, we argue that higher order terms do 
not lead to violation of any of these criteria, and hence that they do not qualitatively change any of our conclusions. 
In particular, the relations between the scaling exponents, and hence the scaling forms ( p7| ) and are still 

correct. The only impact higher order terms have, is to modify the expressions in Eq. for the rcnormalized gaps, 
thereby slightly changing the average ratio a in Eq. ( p3|) and the scaling exponent a. The low-temperature forms of 
the thermodynamic quantities derived in Sec. ^ are valid even though the actual value of the exponent a may shift 
slightly. In particular, it is important to note that the Curie-like form of the magnetic susceptibility does not involve 
a. 



A. Higiier order contributions 

We consider effective interactions between spins which are separated by d (effective) lattice spacings. In the RG 
transformation these are generated by second (and higher) order terms in e, and physically they represent interactions 
mediated by excitations within the locked effective spins. These terms will hence appear even if they are absent in 
the original Hamiltonian. However, the interactions generated in this way fall off exponentially with distance d as 

^-e' . (42) 

as can be seen from the following argument. For given energy scale Aq, consider the strongest bond with spins Sl and 
Sr which are coupled via long-range interactions to the spins Sa and Sb, respectively, with the gaps A^^/Aq ~ e''-*^ 
and Arb/Aq ~ e"*"^ (see Fig. 7). The induced interaction between Sa and Sb due to second-order terms is 
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'^AB _ AalAbB ^ ^dAL+duB _ (43) 

Therefore the non-nearest neighbor interactions decay exponentiaUy with e « 0.2 at the fixed point and cannot lead 
to frustration effects. Consequently, it is justified to restrict our consideration to the dominating nearest-neighbor 
interactions only. The only remaining isotropic interactions are higher order nearest neighbor spin terms (S^ • Si+i)™. 
These terms are local and could in principle be included when calculating the lowest energy spin multiplet and the 
gap to the first excited state in a link. These higher power spin terms are also higher order in e and are unlikely to 
become large enough to change the low-energy spectra of the strongest link qualitatively. 



B. Three-spin decimation 

For some particular combinations of spins and couplings, the renormalized gap becomes larger than the gap just 
removed, A > Ap.g In this section we argue that even in these cases it is justified to use the RG transformation 



outlined in Sec. II B 



For A > Aq a more correct procedure would be to solve the three-spin problem involving the two spins on the 
strongest link and the spin on the link with A > Ao. We would represent the ground state multiplet of the three-spin 
system with one effective spin S and finally calculate the effective couplings between S and its neighbors (Fig. ||b). 
Here we claim that we can obtain essentially the same result using our RSRG scheme (Fig. |^a). In the first step 
the strongest link is replaced by an effective spin S', and the gaps Ai and A2 are renormalized. The renormalized 
gap A2 > Aq by assumption immediately becomes the largest gap in the chain so that, in the next step, the link 
{A2, S", 5*2} is replaced by an effective spin of size S ~ |S" ± 5*21 = \Sl ± Sr ± 521. In this process the gaps A3 
and A']^ are renormalized. The size of the effective spin S in Fig. ^ is given by the absolute value of the (vector) 
sum of the spins S^,, S^j, and S2 parallel or antiparallel according to the sign of the couplings. This is the same 
spin we expect for the ground state multiplet of the three-spin system, i.e., S in Fig. |^b. Similarly we get the sign 
of the couplings between S and its neighbors in Fig. from (^) by aligning the spins according to the signs of the 
couplings and by comparing the direction of the spins Si and S3 with the direction of the effective spin S. The signs 
of the couplings obtained in this way agree with what one expects for the signs of the corresponding couplings in 
the three-spin treatment in Fig. Therefore there can only be a difference in the renormalized coupling strengths 
between the twice- two- spin and three-spin decimation scheme. In both ways, however, the effective couplings Ai and 
A3 are proportional to Ai and A3, respectively (and independent of the magnitude of the unphysically large gap A2 
in the RG treatment). Therefore the discrepancy is minor and will not cause any qualitative difference. 



VII. CONCLUSIONS AND SUMMARY 



We have studied spin chains with random couplings and random spin sizes by means of a real-space RG scheme 
that successively replaces strongly correlated spin-pairs by effective spins. The RG transformation preserves the 
functional form of the Hamiltonian but changes the probability distribution of the links (couplings and spins). This 
procedure generates interactions among the remaining spins. For low enough energies the probability distribution 
of links acquire a scaling form, and for not too singular initial distributions of gaps (so that Ay'=P(A) is regular, 
Uc ~ 0.7) the fixed-point distribution is universal. From a random walk picture for the formation of the large effective 
spins, we argue for a relation between the scaling exponents of length and the average spin size, which together with 
other considerations reduces the number of independent scaling exponents to one. This is confirmed in numerical 
simulations of random spin chains, and in the universal regime, we numerically determine the remaining independent 
scaling exponent to a = 0.22±0.01. 

At low energies (low temperatures) the random spin chain is characterized by large effective spins which interact 
weakly with their nearest neighbors. As temperature is further lowered, the average size of the effective spins increases 
as T~" while the average distance between two effective spins (in units of the original lattice constant) increases as 
rp-2a fpj^jg regime is also characterized by universal temperature dependence of thermodynamic quantities. 

The slow approach to the fixed point (our numerical simulations indicate a crossover region of more than five orders 
of magnitude for reasonable starting configurations) suggests that the true scaling regime may be hard to reach in 
experiments. However, the formation of large effective spins occurs at considerably higher energy scale, and even 
if the scaling exponent a may not have stabilized to its fixed-point value, the distribution of links is roughly like 
the fixed-point distribution. The clearest signal of the formation of large effective spins is perhaps the Curie-like 
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temperature dependence of the uniform magnetic susceptibility, x(r) cx 1/T, in a temperature regime ksT < J 
(J being the typical exchange interaction in the initial spin chain). Also, since the 1/T dependence in x emerges 
before the distribution of links approaches the fixed point, the Curie-like susceptibility should be easier to address 
experimentally than other thermodynamic quantities which may develop scaling behavior only at inaccessibly low 
temperatures. 

When the scaling regime is realized, the most straightforward way to measure the exponent a, and hence the 
rate at which spin degrees of freedom freeze out is through the specific heat, C{T) (x T^"|lnT|. An alternative 
approach which avoids the difficulties connected with measuring small heat transfers at low temperatures is to lower 
the temperature in the presence of a magnetic field until the magnetization saturates. The scaling exponent a may 
then be deduced from the predicted field dependence of the saturated magnetization, Msat{H) cx 7Jt+^. 

The low-energy physics of the random spin chain studied in this paper is very different from that of uniform 
spin systems as well as spin-i random bond anliferromagnetic chains studied in Refs. p"^- pT|Jl9|j20| . In the RAF the 
ground state is a random singlet phase (RSP)E3 where each spin forms a singlet with another spin which may be 
located far away. In the RSP the coupling between two spins that have survived down to some energy scale Aq is 
mediated by virtually exciting all intermediate singlets, leading to effective couplings that decreases exponentially 
with length, J cx exp(~y^). By inverting this relation it follows that length scales logarithmically with the energy, 
n cx I In Aop. For the RSP, arguments analogous to those in Sec. lead to entropy cr(T, H = 0) / L oc \ InTj^^ and 
magnetic susceptibility x{T) / L ix T-'^\lnT\-^ Hln our terminology this corresponds to a — S — wp to logarithmic 
corrections. This is consistent with the interpretation of a in terms of the average renormalization factor a, Eq. (|2^). 
Indeed, in the case of the RAF the perturbative parameter e = Aj/Aq and hence also a goes to zeroE3 near the fixed 
point, implying a = 0. Clearly the RSP is distinct from the large-spin phase (LSP). 

As pointed out in Ref. ^9|, a third possible state of random spin chains at low temperatures is the random dimer 
solid (RDS), which is easily understood within the MDH RG picture. Assuming that e.g., odd links are on average 
slightly stronger than even links in an antiferromagnetic spin-^ chain, the even links are correspondingly more likely 
to be removed. Since the removal of an even (odd) link leaves a renormalized odd (even) link behind, odd links are 
on average renormalized more frequently than even ones, and the separation in energy scale between even and odd 
links becomes more and more pronounced until all the singlets are on even links, the random dimer solid. Unlike the 
RSP and the LSP, thermodynamic quantities in the RDS show non-universal temperature dependence.Ej 

The existence of various low-temperature fixed-points raises the question oLstability of the various phases. As is 
expected from the discussion above, the RSP is unstable towards dimerizationO In contrast, the large-spin fixed point 
is stable towards dimerization. Unlike the RSP in which spins are always removed in pairs, the RG transformation 
in general removes only one spin (i.e., replaces two spins with one effective spin). Hence odd links are turned into 
even on es and vice versa so that dimerization is irrelevant at the large-spin fixed point. From the discussion in 
Sec. Ill it is also clear that both the RSP and the RDS, which are singlet ground states, are unstable towards a 
small fraction of randomly distributed ferromagnetic bonds and/or large spins {S > 1/2). In both cases we expect 
the spin of the ground state to scale with length as S*^ ^ L consistent with the fixed point studied in this paper. 
This has been confirmed in numerical simulations of random antiferromagnetic chains close to the random-singlet 
fixed point with 5% ferromagnetic bonds (chain E in Tab. |) or with 4% of the spins 5=1 (chain G). In all cases 
studied the chains first approach the fixed-point distribution of the RSP by forming singlets through the removal 
of S = 1/2 spins and AF links. However, as the density of higher spins and/or FM links increases, larger effective 
spins start to form and the distribution of links crosses over to the fixed-point distribution of the LSP. Which fixed 
point (either universal or non-universal) is eventually approached depends crucially on the initial distribution. This 
is because the singlet formation is very efficient in decreasing the effective couplings and quickly builds up a singular 
distribution of gaps. Hence, as long as only singlets are formed the distribution rapidly approaches the very singular 
random singlet fixed-point distribution. The degree of the singularity in the gap-distribution at the point where the 
density of 5* 7^ 1/2 and ferromagnetic couplings becomes substantial, determines the behavior of the spin chain. If 
the singularity generated at this point is less than P(A) ~ A~^'= the chain will flow to the universal large-spin fixed 
point, while for stronger singularity the chain approaches one of the non-universal large-spin fixed points. This opens 
up an interesting possibility to access the non-universal fixed points in experiments by starting with a RAF with a 
properly chosen small fraction of FM bonds. 

An interesting open question is what happens away from the Heisenberg point. Fisher extended the MDH RG to 
anisotropic antiferromagnetic spin-i chains to show that also in theiY-regime {Jf < J'f — Jf) they fiow to a random 
singlet fixed point, at least for broad enough initial distributions. til The inclusion of anisotropy in the generalized 
MDH RG and its impact on the large-spin fixed point are interesting problems which we leave for future research. 
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APPENDIX A: DERIVATION OF THE EFFECTIVE HAMILTONIAN 



In this appendix we give a brief derivation of the effective couplings, Eqs. (7) and (9). 



Consider the four-spin Hamiltonian 



where 



1. First-order perturbation theory 



n^Ha + ni (Al) 



Ho = JoSl ■ Sr (A2a) 
H/ = JiSi-Sl + JsSfl-Sa . (A2b) 

We treat Hi as a perturbation to Hq. In the space of the degenerate ground states of Hq, the spins and Sr form 
a state of maximum (minimum) total spin S for Jq < (Jo > 0) while the spins Si and S2 can point in any direction. 
The degenerate ground states span the Hilbert space H, the product space of the spin spaces for Si, S and S2. Each 
state \m1Mm2) = \mi) (8) \M) (g) |m2) in H is labeled by the corresponding azimuthal quantum numbers mi, M and 
1712- Ti-i partly liftSpihe degeneracy and induces an effective Hamiltonian H'^^ in H. To order J1.2/J0 the matrix 
elements of Ti"^ areN 

^miM™„™;M'™i = {miMm2\ni\m[M'm'2) . (A3) 
We calculate in two steps: We first establish the operator identities (valid in H) 

51- Sl^c{Sl,Sr,S)Si-S (A4a) 

52- Sr = c{Sr,Sl,S)S2-S (A4b) 



which, together with Eqs. (A2b) and (|A3|), give 

= JiSi • S + J2S • S2 (A5) 

with 

Ji^ .hc{SL,SR,S) (A6a) 
J2 = Jic{Sr,Sl,S) . (A6b) 



We then determine the constant c{Sl, Sr, S) by calculating a suitable matrix element of the two operators in Eq. ( [A4a| ) . 
To establish ( [A4a| ) we define the usual raising and lowering operators which act on a state \M) in a spin-S" multiplet 



in the following way. 



^-IM + I) =Ali\M) 
If. I 



s+\m)^aIj\m- 



where the constants Afj = y/ {S — M){S + A/ + 1) assure that (M|M) = 1. Consider first the operators S and 



5j Sj^ : 
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{miMm2\S+Sl\m[M'm'^) = 5m'^,rm-ihi' M+i5m'^,m,A^r^,-i{M\Sl\M + 1) . 
The matrix element in the right hand side of Eq. ( A7b ) is 

{M\Sl\M + 1) = (^f,+i^f,+2 • ■■Al_,)-^M\Sl{S- f-^'-'\S) 

^AlMl-irHs~i\si\s) , 

where we have used the fact that S~ = + commutes with S^. From Eqs. (A7a), (A7b) and (A8) we see that 

S+S^ - {Al_,)-^S-l\SZ\S)S+S- (A9) 

for every combination of mi...M', i.e., as an operator identity in H. We note that (5*— is a real number, and 

take conjugation of Eq. ([A9|): 



(A7a) 
(A7b) 



(A8) 



S^S+ = {Al_,)-\S-l\SZ\S)S^S^ 



We obtain the operator identity 



s^s'i + Si Si = {A%^^)-'{s ^i\si\s)[sis'' + sisy 



(AlO) 



Equation (A4a) then follows from Eq. (AlO) and rotational invariance. From the relation {S — l\Sj^\S) — 
{Al_^)~HS\S+Sl\S) = 2(Af_J-i(5|5£|5) we obtain 



c{Sl,Sr, S) 



{S\Sl\S) 



(All) 



To determine {S\Sl\S) in Eq. ( All ) we evaluate the matrix element (S'ISl • SjS*) in two different ways. First we use 
S = Sl + Sij to get 

{S\Sl-^\S) ^{S\SI\S) + {S\Sl-^r\S) 



- [S{S + 1) + Sl{Sl + 1) - Sr{Sr + 1)] 



The same matrix element can also be written as 



{S\Sl ■ S\S) = {S\SIS^\S) + \{S\SIS+\S) + ]^{s\sts-\s) 



(A12) 



(A13) 



The first term in Eq. ( A13 ) equals S'(S'|S'£|S'), the second term vanishes since \S) is a highest weight state, and for 
the same reason we can replace the operator S^S~ in the third term by the commutator [5^, S^] = 2S'£ so that 



{S\Sl-S\S) = {S + 1){S\SI\S) 



From Eqs. ( |A12[ ), ( |A14D , and (|AllD finally follows 

SiS + 1) + SLiSL + 1) - SniSR + 1) 



c{Sl, Sr, S) 



25(5 + 1) 



(A14) 
(A15) 



For the three cases of interest to us, Eqs. ( A15 ) and ( A6) give 

Sl 



Ji 
Ji 



and 



J2 



Sl + Sr 
Sl + 1 
Sl-Sr + 1 
J 'S'l 
Sr — Sl + ^ 

Sr 



J2 = <^ 



Sl + Sr 

J Sr, 



J2 



Sl — Sr + 1 
Sr + 1 



Jo <0 

Jo > 0, 5l > Sr 

Jo > 0, 5l < Sr 

Jo <0 

Jo > 0, Sl > Sr 

Jo > 0, 5l < Sr 



(A16) 



(A17) 



' Sr — 5l + 1 

from which Eqs. (0) follow by using the relation between couplings J and gaps A, Eqs. 
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2. The effective coupling in the case of singlet formation: generalization of the MDH RG transformation 



In the case where Jq > and Sl = Sr, the effective spin is a singlet, S = 0, and the effective Hamihonian, Eq. 
vanishes. To get a nonzero couphng between and 5*2, we have to include second order perturbation. In the case 
Sl = Sr = ^ this gives the effective couplings_Lised in the MDH RG. Since the first order contribution is zero, we 
have to solve the second-order secular equationcS 

Det (Vgg, - E<^^1Sgg,) = (Al8) 



with 

where \g) is a ground state of Hq, the sum is over all excited states |e) and Eg is the (unperturbed) energy of the 
state |e). Since there is no coupling between the outer spins if either Ji or J2 is zero, only terms proportional to J1J2 
in Eq. ( A19| ) contribute to the effective coupling (the terms proportional to and J| give an overall energy shift 



which, for our purposes, can be discarded). Denoting the states formed by the spins Si and 5*2 by \A), the states 
formed by Sn and Sl by \B) and in particular the ground state singlet by \B = 0) we have (in obvious notation) 



Vaa = 2Ji J2 



(A, 0\SlSi\A', B') {A\ B'\SliSI\A, 0) 
^ -Jo^Sb'{Sb' + 1) 



A'M' ""^-i 

J1J2 V- (A, 0|5^5+|A^ B'){A',B'\S^S+\A, 0) 

2 -■hlSB'iSB' + l) 

J1J2 V- {A,0\S+SZ\A',B'){A',B'\S+S^\A,0) 
2 j^, -Jo^Sb'{Sb' + 1) 

The sum over the states \A') can be performed separately to give 

W2,,,...,.,,^ {0\Sl\B'){B'\S^j,\0) 
= -- \f SB'iSB^ + l) 

'^'•^\A\StS-\A)^ ^^\Sl\B'){B'\S+\Q) 



(A20) 



Jo ' ' ' V Sb-{Sb' + 1) 



The matrix elements in the sums get contributions only from spin-1 multiplets, and hence we can extract the factor 
^/[Sb'{Sb' + 1)] = 1/2 and perform the sum over the complete set of states \B'). Finally, since 

miSrn = ^(0|5+^^|0) = \{Q\S-J%\Q) = i(0|Si • S«|0) = -\sl{Sl + 1) (A22) 

we have 

Vaa = ''JlMl^^p^{A\S, ■ SM) (A23) 

from which it is clear that the effective Hamiltonian in this case is 

U"^ = JSi • S2 (A24) 

with 

j^ 2J,J2Sl{Sl + 1) ^^25) 
3Jo 

In the special case of S'l = S'/j = 1/2 this gives J = Ji J2/2J0 in agreement with Ref. Ilj. 
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APPENDIX B: THE GENERALIZATION OF THE MDH RG FLOW EQUATIONS 



Although in this paper we have studied the RG flow by numerical simulations, it is of interest to show that flow 
equations of the form Eqs. (^3|) exist and can be expected to have attractive fixed point(s) to which the probability 
distributions eventually flow. Rather than deriving the explicit forms of Fi and i^2, which are rather complicated, we 
will outline how this can be done and point out to what extent the generalized flow equations ( [T3| ) differ from the 
MDH RG flow equation (|). 

We may view one step in the RG transformation as the removal of one link {Aq, Sl, Sr} and the change of two 
links, {Ai, Si,Sl}, {A2, Sr, S2} — > {Ai, S*!, S}, {A2, S, 82}- For an infinitely long chain we remove a small fraction 
of links with gaps in the interval [Aq — dAo, Aq]. This is illustrated in Fig. |l^, where, for convenience, we show 
P"*(Ao; A, Sl, Sji) and P^{Ao; —A, Sl, Sr) in the same diagram. At an energy scale Aq, and are non-zero 
only for links with A < Aq. We remove all links in the thin shell A e [Aq — dAo, Aq], which causes a small fraction 
of links inside the equal-gap surface A = Aq to hop around. The changes in P^ and P^ due to the links that move 
around are of order dAo so that in the limit dAo — > the RG flow is described by a set of first-order differential 
equations 

Fi[P^,P% (Bla) 

F2[P^,P% (Bib) 

where Fi and F2 are two (non-linear) functionals of P^ and P^ , whose explicit forms depend on the functions /„ in 
Eq. (^). If renormalized gaps were always smaller than Aq, it would be straightforward to write down the explicit form 
of the flow equations (O). Assuming no correlations between neighboring links (except for the obvious correlation 
that they share one spinj^ we find four types of terms in Fi; (i) One term proportional to —S{A — Ao)P'^ that depletes 
the region of links where A £ [Aq — dAo, Aq]. (**) One term that decreases P^ due to links {A,S'i,S'l} that are 
transformed because they neighbor a link that is replaced by an effective spin, (in) A set of terms that increase 
P^{Ao; Ai, S*!, S") because some links are transformed into links {Ai, ^i, 5}. (iv) One term proportional to P^ that 
compensates for the overall decrease in the number of links and keeps the probability distributions normalized. The 
functional F2 has a similar structure. Only the third kind of terms involve the RG functions /„. We note that in the 
original MDH RG the terms (ii) and (iv) cancel. The term (i) can be omitted if the probability distribution is defined 
only on the interval [0, Aq]. This was done in Ref. |lj so that the only term that appears in Eq. (^) is one term of 
type [iii). 

If some renormalized gaps become larger than Aq, then before we take the limit dAo — > we must consider the 
fraction of the transformed links that acquire A > Aq. As we have discussed in Sec. VIB, in our RG scheme a finite 
fraction of the links acquires a larger gap A > Aq; in Fig. |l^ this corresponds to links that jump outside the support 
of P^ and P^ . If the unphysically strong links are not taken care of before the next slice of links are removed, more 
and more links will end up outside the Aq surface. In the limit dAo this will make it impossible to define Aq, 
since smoothly integrating out links in a finite gap interval in this case generates a small but finite probability for 
gaps of any strength. This problem is remedied in the following way. Since in our RG scheme, no property of the 
renormalized links depends on the actual value of the gap in the strongest link, we can modify the equations (^ so 
that A = min{A/„, Aq}. This does not change the discrete RG scheme since the renormalized link is removed in the 
next step of the RG after which there is no information in the chain of the value of A. In Fig. |l^, this modification 
of the functions /„ is illustrated by projecting strong links back onto the surface A — Aq. The projected links are 
then removed, which transforms some links in P^ and P^ , and also generates an even smaller fraction of links above 
the Ag-surface. These links are then projected and removed etc. Keeping terms of order dAo, the infinite sequence 
of removing smaller and smaller fractions of links that in the previous step have jumped out of the distribution 
contributes a series of smaller and smaller terms in the shift of P'^ and P^ that has to be summed before the limit 
dAo — > is taken. This resummation makes the functionals Fi_2 in ( p^ ) rather complicated. 



dP^ 

dA^ 

dpP 

'dA^ 
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TABLE I. The initial conditions for the 16 chains simulated numerically. 



FIG. 1. Schematic pictures of the RG scheme, (a) The original MDH decimation, (b) Definition of a link as two neighboring 
spins Sl and Sr and the gap A. (c) The generalized MDH decimation. 

FIG. 2. (a) The antiferromagnetic fixed-point distribution of spins and gaps, Qsa(A/Ao,5) = JJ^ d5'Q'^(A/Ao, 5, 5'). 
The spins are in units of (5) and the distribution is normalized according to (hll) . (b) The ferromagnetic fixed-point distribution 
of spins and gaps, Qg/^{A/ Aq, 5), defined analogously to Qsa in {^)- (c) The antiferromagnetic fixed-point distributions of left 
and right spins, Q'§g{S, S') = dxQ^{x, S, S'). The units and normalization are as in (a), (d) The ferromagnetic fixed-point 
distributions of left and right spins, Qgg{S, S'), defined analogously to Qgs in (c) ^nd with units and normalization as in (a). 

FIG. 3. (a) The distributions of gaps at the fixed point, Q^'^(A/Ao) = dSdS'Q^-''{A/Ao, 5, S'). (b) The distributions 
of spins at the fixed point, Qg'^(5) = J^^ Jq° dS'Q^'^ {x, 5, 5'). The spins are in units of (5). 

FIG. 4. (a) The average gap (A) as a function of Ao in chain C in Tab. |. (b) The average spin (5) and length n as a 
function of Ao in chain C in Tab. Ill 



FIG. 5. Efi'ective exponents as functions of Aq ^ in chain C in Tab. 



FIG. 6. The fixed-point distributions of AF gaps for four singular chains (chains E, G, K and O in Tab. |) (dotted lines). 
These are to be compared with the corresponding fixed-point distribution of regular chains (solid line) 



FIG. 7. Log-log plot of the distribution of gaps at the fixed point for regular initial distributions. 



FIG. 8. Higher order terms induce a coupling between Sa and Sb as Sl and Sr are replaced by an efi^ective spin. 

FIG. 9. (a) The two-step decimation of three spins used in the numerical simulations, (b) The three-spin decimation in one 
step. 
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FIG. 10. The link distributions are non-zero only for A < Aq. Replacing gaps in the thin slice A £ [Aq — dAg, Aq] shifts a 
fraction of links in and . 
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Figure 2a in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 




Figure 2b in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 




Figure 2c in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 
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Figure 2d in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 
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Figure 3 in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 




Figure 4 in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 




Figure 5 in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 
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Figure 6 in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 




Figure 7 in the article "The low-energy fixed points of random 
quantum spin chains" by E. Westerberg, A. Furusaki, M. Sigrist, 
and P. A. Lee. 
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Figure 10 in the article "The Low-Energy Fixed Points of Random Quantum 
Spin Chains" by E. Westerberg, A. Furusaki, M. Sigrist, and P. A. Lee 




